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ABSTRACT 

We present a catalog of compact sources derived from the QUaD Galactic Plane Survey. The survey covers 
^ 800 square degrees of the inner galaxy {\b\ < 4°) in Stokes /, Q, and U parameters at 100 and 150 GHz, 
with angular resolution 5 and 3.5 arcminutes respectively. 505 unique sources are identified in /, of which 239 
are spatially matched between frequency bands, with 50 (216) detected at 100 (150) GHz alone; 182 sources 
are identified as ultracompact H II (UCH ll) regions. Approximating the distribution of total intensity source 
fluxes as a power-law, we find a slope of 75,100 = -1 .8 ± 0.4 at 100 GHz, and 75,150 = -2.2 ± 0.4 at 150 GHz. 
Similarly, the power-law index of the source two-point angular correlation function is 7^,100 = -1 .21 ± 0.04 and 
79.150 = -1.25 ±0.04. The total intensity spectral index distribution peaks at a/ ^ 0.25, indicating that dust 
emission is not the only source of radiation produced by these objects between 100 and 150 GHz; free-free 
radiation is likely significant in the 100 GHz band. Four sources are detected in polarized intensity P, of which 
three have matching counterparts in /. Three of the polarized sources lie close to the galactic center, Sagittarius 
A*, Sagittarius B2 and the Galactic Radio Arc, while the fourth is RCW 49, a bright H II region. An extended 
polarized source, undetected by the source extraction algorithm on account of its ^ 0.5° size, is identified 
visually, and is an isolated example of large-scale polarized emission oriented distinctly from the bulk galactic 
dust polarization. 

Subject headings: Surveys — Galaxies: Milky Way — Galaxy: Structure — ISM: star formation — H II 
regions 



1. INTRODUCTION 

Millimeter (mm), sub-millimeter (sub-mm) and far-infrared 
(FIR) observations are ideal for studying the properties of 
star-forming regions in the galaxy, in particular the cool en- 
velopes of dust and gas which host sites of potential and active 
star formation. By spanning the peak in the spectra of these 
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objects, measurements between the mm and FIR can tightly 
constrain the parameters of the thermal radiation produced by 
the dust. In particular, the mass of a star-forming core and 
its surrounding envelope is well-traced by its measured flux 
in these bands, since this radiation is optically thin at sub-mm 
and longer wavelengths. 

Surveys covering large sections of the galaxy have the po- 
tential to collect statistical samples of cores in a range of evo- 
lutionary states, comparatively free of bias introduced by tar- 
getting particular regions. These surveys are ideal to study 
processes related to star-forming regions, such as measuring 
the core mass function (from which the stellar initial mass 
function may be derived), particularly at the high-mass end, 
which, on account of the short-lived high mass cores, is un- 
derstudied relative to lower masses (e.g. [Enoch et al.|2006[ 
Young et al.'2006; Enoch et al. 2008). Combination with in- 



frared (IR) data yields insight into the ages of cores, permit- 
ting differentiation between prestellar sub-mm cores, which 
lack an IR counterpart, and protostellar cores, in which the 
ultraviolet radiation produced by protostars is re-radiated into 
the mm, sub-mm and IR by the surrounding envelope. Phe- 
nomena associated with later evolutionary phases, such as 
mass ejection, dissipation of the envelope, and dynamical in- 
teractions are not significant in the prestellar or protostellar 
stage — the mass and spatial distribution of such cores there- 
fore capture information regarding the fragmentation pro- 
cess (Enoch et al. 2006). 

Observations of polarized radiation permit a window to 
study the role of magnetic fields (e.g. [Greaves et al.||T995| 
Novak et al. 1997| l, and their role in providing support against 



collapse. In the mm and sub-mm, polarization is due to emis- 
sion along the long axis of dust grains partially aligned by 
the magnetic field, and thus measurements of the dust polar- 
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Fig. 1 . — Fourth (left) and third (right) quadrant field-differenced 100 GHz Stokes / map smoothed to the beam scale (5'), with color scale in MJy/sr. The solid 
black lines indicate the survey coverage. 



ization directly probe local magnetic fields (e.g. |Hildebrand 
[1988) . These fields are thought to strongly influence the evo- 
lution of molecular clouds, since they provide support pre- 
venting the collapse of the gas and subsequent triggering of 
star formation. 

Several large-scale surveys are akeady underway or com- 
pleted to help address these questions. Herschel ( Pilbratt] 
etal.|20 10) and Planck ( |ViUa et al.|2002| l will provide exten- 



sive spectral coverage from the radio to the far infrared, fully 
characterizing the spectral energy distribution (SED) of star- 
forming cores over the fully sky; selected existing results in 
targetted regions include (e.g. Hennemann et al.|2010 Juvela 



|et al. 2010 Andre et al.|2010| l, but are limited to total intensity 
observations. Ground and balloon instruments also contribute 
substantially to the literature: SchuUer et al. ( 2009 [ l present 
an APEX LABOCA 95 deg^ survey in total intensity with res- 
olution of 19.2" at 353 GHz, with the final survey coverage 
expected to reach 350 deg^; Bolocam has mapped 150 deg^ 
of the first galactic quadrant at 1.1mm (268 GHz) with resolu- 
tion 33", with a source catalog pr esented in|Rosolowsky et"ar 
( [20T0| ); BLAST ( jOlmi et al.|2()09t |Netterfield et al.|2009 ') pro^ 
vide a 50 deg^ survey of the Vela molecular cloud at 250, 
350 and 500 microns (36, 42 and 60 arcsec resolution respec- 
tively). 

Observations at comparable resolution are currently scarce 
at ^100 GHz, and yet provide additional constraining power 
to the Rayleigh- Jeans tail of the thermal dust spectrum, and 
probe for contributions due to other emission mechanisms 
which contribute increasingly at lower frequencies (e.g. free- 
free). Furthermore, there is little high angular resolution po- 
larization data at these frequencies, despite their utility in un- 
derstanding star-forming regions. 

In this paper we present a catalog of compact sources found 
in the QUaD galactic plane survey ( jCulverhouse et al.|2010| l, 
which covers over ~ 800 square degrees of the low-latitude 
galactic plane at 100 and 150 GHz with beam FWHM of 5' 
and 3.5' respectively, in Stokes /, Q and U parameters A 
survey of this size, frequency and angular resolution can be 
used to investigate the polarized and unpolarized properties 

' The QUaD maps and source catalogs analyzed 
in this paper are available for public download at 

http :/ /find . spa . umn . edu : /quad/quad_galactic/ 



of both diffuse emission and discrete sources. The QUaD sur- 
vey was conducted blind, in that no region was specifically 
targetted. In principle, this allows the construction of statisti- 
cal samples of cores, representative of the distribution of core 
masses and ages in the galaxy as a whole. However, we note 
that at the ^few arcminute resolution of the survey, the maps 
do not generally resolve individual cores: Dense cores typi- 
cally have size ^ O.lpc (e.g. Williams et al. 2000), hence for 
nominal distances of a few hundred pc, the sources presented 
here should be considered as 'clumps' hosting cores rather 
than individual cores themselves. In addition to the lack of 
resolution and accurate clump distances, the contribution of 
free-free emission at 100 GHz biases measurements of the 
flux from the dust component; these caveats prevent reliable 
mass calculation. Our goals here are therefore to analyze the 
observed quantities of the sources in the survey, rather than 
infer their physical properties. 

Basic information on the instrument, observations and maps 
is presented in Section |2] In Section [3] we describe our algo- 
rithm for extracting sources in the presence of a diffuse back- 
ground. The global properties of the catalog are discussed in 
Section |4] with the full catalogs presented in Table |3] (total 
intensity) and Table |4] (polarized intensity). Discussion and 
Conclusions are found in Section |5] Extensive simulations, 
presented in Appendix [A] are used to quantify the effects of 
mapmaking and source extraction algorithm on the recovered 
source properties. 

2. INSTRUMENT, OBSERVATIONS AND MAPS 

Here we summarize the features of the QUaD galactic plane 
survey. A detailed description of the instrument can be found 
in |Hinderks et al.| ( [2009) , hereafter referred to as the "Instru- 
ment Paper". The field selection, survey strategy, data pro- 
cessing and construction of the Stokes /, Q and U maps are 
presented in Culverhouse et al. (2010 1, hereafter the "Map 
Paper". 

QUaD was a 2.6 m Cassegrain radio telescope on the mount 
originally constructed for the DASI experiment (Leitch et al. 
|2002( l. This is an az/el mount, with a third axis allowing the 
entire optics and receiver to be rotated around the line of sight. 
The mount is enclosed in a reflective ground shield, extended 
from DASI, on top of a tower at the MAPO observatory ap- 
proximately 1 km from the geographic South Pole. 



QUaD Galactic Plane - Compact Sources 



3 




model map 







-55.5' 



240 



238 
RA (deg) 



236 



-55.5 



240 



238 
RA (deg) 



236 



55.5' 



240 



238 
RA (deg) 



236 



-55.5' 



240 



-53.5 



^ -54 



^ -54.5 



-55.5 




238 
RA (deg) 

m - model 



236 



240 238 

RA (deg) 



-53.5 



?P -54 



^ -54.5 



-55.5 




240 238 

RA (deg) 



240 238 236 

RA (deg) 



-53.5 



?P -54 



^ -54.5 



-55.5 




240 238 236 

RA (deg) 



Fig. 2. — Source segmentation and fitting using a section of QUaD 100 GHz / data; colorscale is the same in all panels and runs from -0.0035 to 0.035 MJy/sr. 
Top Left: Raw map nir. Bottom Left: Initial estimate of background map nii^j. Top row, second column: Initial background subtracted map m,. Bottom row, 
second column: Background map corrected for discrete source flux (miyi-i^j-orr)- Top row, third column: Background subtracted map m, using source-corrected 
background map nitck.corr- Bottom row, third column: Map used for source fitting mr—mi,ck.corr- Top right: Model of discrete source population using fits to mf. 
Bottom right: Residual between model sky and input image nir. 

The QUaD receiver consisted of 31 pairs of polarization 
sensitive bolometers (PSBs; |Jones et al.|2 003 ), 12 at 100 GHz, 
and 19 at 150 GHz. The bolometers were read out using AC 
bias electronics, and digitized by a 100 Hz, 16 bit ADC; the 
raw data were staged on disk at Pole and transferred out daily 
via satellite. 

The observations reported in this paper were made between 
July and October 2007, with the telescope decommissioned 
in late 2007. In total QUaD surveyed the galaxy for 40 days, 
covermg a total of - 800 deg^. The survey is divided into 
two regions, approximately covering 245-295° and 315-5° in 
galactic longitude Z, and -4 to H-4° in galactic latitude b. These 
regions are loosely termed the 'third quadrant' and 'fourth 
quadrant' throughout. 

Maps are made by coadding the timestream from each de- 
tector into flat-sky (RA, dec) pixels of size 0.02° x 0.02°. 
The absolute pointing accuracy was determined to be ^ 0.5' 
rms, using pointing checks on RCW 38 and other galactic 
sources over two seasons of CMB observations (see Instru- 
ment Paper for further details). A field-differencing scheme 
was used to remove spurious ground contamination; all the 
results presented here are derived using field-differenced 
maps. Absolute calibration is applied using a scaling fac- 
tor at each frequency, derived in the QUaD CMB analysis 
presented in [Brown et al.| ( |2009 1 



Stokes parameter are also constructed, which give a measure 
of the noise in each map pixel. The typical survey sensitiv- 
ity in each survey area is 74 (107) kJy/sr at 100 (150) GHz 
in I, and 98 (120) kJy/sr in QIU, at a spatial resolution of 5 
(3.5) arcminutes at 100 (150) GHz. The orientation of Q and 
U in the QUaD polarization maps follows the lAU conven- 
tion ( [Hamaker & Bre gman 1996), in which +Q is parallel to 
N-S and +U parallel to NE-SW. 

3. SOURCE EXTRACTION 

As is readily apparent from Figure [T] a substantial contrib- 
utor to the sky signal is diffuse emission. This 'background' 
increases the uncertainty in measured properties of compact 
sources above that due to detector and atmospheric noise. 
However, the systematic effect of diffuse emission can be re- 
duced using spatial filtering. 

The source extraction method implemented here is an adap- 
tation of the algorithm described in Desert et al. ( 2008| l (a 
mexican-hat linear filter in image space), which was designed 
to separate compact sources from the diffuse galactic emis- 
sion for the Archeops experiment. Maps of each stage of 
the source extraction algorithm, described below, are shown 
in Figure [2] in a representative section of the fourth quadrant 
100 GHz /data. 



These factors were calcu- 
lated by cross-calibrating QUaD CMB maps to those from 
the Boomerang experiment (Masi et al. 2006 1, and have an 
estimated uncertainty of 3.5%. The 100 GHz / map for both 
survey regions is shown in Figure [T} the reader is referred to 
the Map Paper for similar maps at both frequencies and in /, 
Q, and U. In addition to the sky maps, variance maps for each 



In our algorithm, two sets of smoothed / and P = \/Q^ + U^ 
maps are made at each frequency; both are derived from the 
raw maps m^. The first set, m^, consists of smoothed to the 
beam scale abeam- In the second, nir is smoothed to an angular 
scale ahck = 2.5 x abeam to form a template map of the diffuse 
background, nibck- In both cases, a circularly symmetric gaus- 
sian function is the smoothing kernel. The choice of abck is 
designed to minimize the background contribution to source 



4 



QUaD collaboration 



fluxes without introducing large biases in the measured flux. 
In Appendix |A.5| simulations demonstrate the consequences 
of other choices of atck on recovered source fluxes. The back- 
ground maps are then subtracted from beam-smoothed maps 
to yield the 'source extraction' map m,: 

m:,, = mi,t-mi.hck (1) 

and likewise for P. Note that since <Jbck > o'beam. Equation [T] 
is equivalent to convolving with filter constructed from the 
difference of the two smoothing kernels, commonly referred 
to as a 'mexican hat' filter. 

Negative pixels due to ringing are masked, and the remain- 
ing pixels subjected to signal-to-noise thresholding. Pixels 
above a signal-to-noise threshold of 5 (3) for total (polarized) 
intensity are flagged as belonging to source candidates. The 
polarization data has a lower extraction threshold because the 
noise properties are closer to white on account of the unpo- 
larized atmosphere, and also because the diffuse component 
amplitude (fractional polarization < 2%; see Map Paper) is 
close to the instrumental noise level and therefore its effect 
on source fluxes is small. High signal-to-noise regions in 
the P map define a set of pixels to which we fit polarized 
sources in the Q and U maps separately. Candidate pixels 
in all Stokes maps are subject to suitability checks; isolated 
pixels or groups of pixels smaller than the beam width are 
removed. 

Sources in close proximity tend to be members of the same 
thresholded region, so an internally-developed segmentation 
algorithm based on the SExtractor code ( jBertin & Arnouts] 
|1996| l is used to split such regions into separate sources. 
Source segmentation is applied to the /, Q, and U maps sepa- 
rately at each frequency, resulting in a set of six source posi- 
tion lists, along with the map pixels assigned to each source. 

Having determined source positions, the background maps 
mhck are regenerated by again smoothing nir, but with the pix- 
els corresponding to discrete sources replaced by the local 
median — this 'source-corrected' background map is denoted 
mck.corr, with t7i,_corr = " mtck,corr following from Equa- 
tion [T] The median-replacement step reduces the amount of 
ringing due to the background filtering implemented in Equa- 
tion[T](see Figure[2]i. The resulting background map contains 
less leaked flux due to smoothing discrete souces with a ker- 
nel larger than the beam size. The background subtraction, 
source detection and segmentation stages are then repeated. 

All ingredients for measuring source properties are present 
at this point: the background map mhck.mrr, a list of pixels be- 
longing to each source, and the input map itself m, , with its 
variance map <j^. The background map nihck.corr is subtracted 
from the input map nir yielding the map to which source mod- 
els are fit, 

mij = mi.,r-m,.hck,corr- (2) 

and likewise for P. An elliptical gaussian is fit to the resulting 
pixels for each source; an example of the model reconstructed 
from these fits is shown in Figure [2l Residuals of this model 
against m^, also shown in Figure [2f indicate that the source- 
fitting works well, except for two cases: 1) particularly bright 
sources, which can leave residuals at the ^ few percent level; 
and 2) sources in close proximity, for which the source seg- 
mentation fails and the sources are classed as a single object. 

Fits are performed independently at each frequency in /, Q 
and U, with pixel noise taken from the corresponding vari- 
ance map a^; errors on source properties follow directly from 
the parameter errors returned by the fit minimizer. The el- 
liptical gaussian fit parameters are used to calculate source 



properties such as flux and position; in a small number of 
cases, the uncertainty on a fit parameter diverges, in which 
case we do not quote the uncertainty on any physical quan- 
tity derived from this parameter Derived quantities such as 
spectral index a = \og{l2/I\)/\og{v2/v\), polarization angle 
(f) = O.Stan"' {U /Q), polarization fraction f//, and their asso- 
ciated errors are calculated from /, Q and U fluxes. In general, 
a source detected in / will not have the same set of pixels as 
in Q or U as each map is treated independently; sources are 
spatially matched across catalogs later to determine e.g. po- 
larization fraction. Given two catalogs A and B (such as total 
intensity at two frequencies), each source in A is matched to 
a source in B if their angular separation is less than three map 
pixels (3.6'), conservatively larger than the rms day-to-day 
telescope pointing wander of ^ 0.5'. If more than one source 
in B matches a source in A, as can happen when matching 
sources between 100 and 150 GHz due to the higher resolu- 
tion in the latter band, the closer of the two is selected. Having 
matched sources, the corresponding physical quantities are 
combined to yield the derived quantity such as the spectral 
index. 

3.1. Consequences of Field-dijferencing 
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Fig. 3. — Top: Subsection of field-differenced 100 GHz / map cover- 
ing galactic center. Few sources lie near the left (high R.A.) edge of the 
map, where the trail field is aligned. Bottom: Subsection of field-differenced 
100 GHz / map covering the lowest Decl. of the survey. Two real sources 
close to the trail field region (high R.A. edge) are circled with solid black 
lines. Their field-differenced counterparts are circled with broken black lines; 
notice the change in sign in intensity of these spurious sources (other field- 
differenced sources are also visibile). No sources lie as close to the trail field 
edge in the top plot, and hence no spurious sources are detected at this Decl. 

The field-differencing operation performed on the 
timestream to remove ground contamination can result in 
spurious sources in the data. If a source lies in the trail 
field of the observations (at larger R.A.), when the trail field 
is differenced against the lead field the source will appear 
negative in the lead field, resulting in a negative measured 
flux. Such sources are removed in / by rejecting candidates 
with fluxes below zero at a signal-to-noise of 5 or greater. 
This rejection is not possible in 2 or f/ because the polarized 
flux can take positive or negative values. Instead, detected 
polarized sources are matched to the total intensity source 
catalog; if the polarized source is matched to a source of 
negative total intensity with |S/N| > 5, the polarized source 
is removed from the catalog. This method may allow small 
numbers of field differenced polarized sources to leak into 
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TABLE 1 
Source Statistics 



Type 


100 GHz 


150 GHz 


Freq Matched 


/ 


289 


455 


239 


P 


3 


3 


2 


IIP Matched 


2 


2 


1 



the catalog, since sources in / are extracted at a higher signif- 
icance threshold {5a) than polarization — a field differenced 
source in P may not be matched to a total intensity candidate 
and therefore cannot be rejected. 

Field-differenced sources are also increasingly expected at 
low declinations. The central PSB pair are aligned on the 
plane of the galaxy, with the low R.A. edge of the trail field 
aligned with the high R.A. edge of the lead field; the width 
of the QUaD focal plane allows for overlapping coverage of 
the lead and trail fields. At higher elevation (lower declina- 
tion) the scan throw of 5Az =15° translates into a smaller R.A. 
range as JR. A. = 6Az x cos(Decl.). Due to the alignment of 
lead and trail field edges and the decreasing scan throw in 
R.A. at lower Decl., the trail field lies closer to the galactic 
plane at low Decl. — see Figure[3](we further refer the reader 
to Figure 2 of the Map Paper for a graphical representation of 
the lead/trail field geometry over the full survey). At lower 
Decl. the trail field is therefore more likely to contain a bright 
source (Figure |3]l. However, the QUaD catalog indicates that 
most detected sources lie within three degrees of the galactic 
plane (Section [4.2| i, and thus the contribution of field differ- 
enced sources is small over most of the survey, since the trail 
fields never get closer than (15°/2) x cos(60) = 3.75° to 
b = Q. 

4. RESULTS 

Source catalogs from both the third and fourth quadrant 
maps are extracted and combined for the purposes of calculat- 
ing statistical properties. Spurious sources are rejected if the 
flux is negative at one frequency and undetected at the other 
frequency. Statistics for the survey are shown in Table [T] 

In total, 289 (455) sources are detected in / at 100 (150) 
GHz, of which 239 are spatially matched between frequency 
bands, resulting in 505 unique sources in total intensity. Po- 
sition, major/minor axes, flux, and spectral index for each 
source in / are given in Table [3] Four sources are detected 
in P, of which two are matched spatially across bands; three 
of these polarized sources have matching counterparts in at 
least one frequency band in /. Properties of these sources are 
presented in Table |4] 

Simulations are used to determine survey completeness and 
purity, and accuracy of recovered source parameters. Four dif- 
ferent types of simulations (labelled Siml-Sim4) are used to 
calculate these quantities and how they are affected by fea- 
tures particular to a galactic plane survey, specifically the 
influence of the diffuse background, and the effect of an 
anisotropic distribution of spatially clustered point sources. 
A detailed description of the simulations is presented in Ap- 
pendix |A] 

The 90% survey completeness in total intensity, C/ go, is de- 
termined from Sim4 (the most realistic simulation used, incor- 
porating both detector, atmospheric noise, spatially clustered 
point sources and a model for the diffuse background); we 
find 5.9 and 2.9 Jy at 100 and 150 GHz respectively. These 
do not change significantly with the inclusion/exclusion of a 



diffuse background, indicating that the background removal 
strategy described in Section[3]is effective. 

In polarization, the 90% completeness Cptjo = 20.3 Jy at 
100 GHz, and 1.1 Jy at 150 GHz. The large difference is 
likely due to confusion when projecting source emission of 
randomly oriented polarization angle along the line of sight, 
an effect reduced at 150 GHz on account of the smaller beam 
size; randomly distributed sources give a 9 0% c onfusion limit 
of ^ 1 Jy at both frequencies (Appendix A. 2 1. Similarly to 



total intensity, the completeness limit in polarization is not 
strongly dependent on the presence of diffuse galactic emis- 
sion. 

At the signal-to-noise extraction threshold of 5 in /, the sur- 
vey is 90- 100% pure at both frequencies, while the purity 
in polarization is very similar, but at an extraction signal-to- 
noise of 3; the spread in purity arises due to the different types 
of simulation. 

The QUaD catalog is matched to the IRAS Point Source 
Catalo g (IRAS-PSC) and the PMN catalog ( [Condon et aT] 
19931 1 using a search radius of 0.1'. Of the sources detected 
at 100 (150) GHz, 97% (87%) have IRAS-PSC counterparts; 
this discrepency could be due to the similar resolution of 
IRAS lOO^m and QUaD 100 GHz, while sources in close 
proximity may be resolved at 150 GHz due to the higher an- 
gular resolution. Three sources without IRAS-PSC counter- 
part have associations with PMN sources. That such a large 
fraction of sources have IR counterparts indicates that cores 
located inside the detected clumps are past the prestellar phase 
and have thermally radiating dust envelopes. This result might 
be expected; since the QUaD frequency bands probe the dust 
emission well away from the core SED peak, we are unlikely 
to detect the prestellar or starless cores which consist solely 
of very cold molecular gas and have no internal source of lu- 
minosity. 

The QUaD catalog can provide constraints on the con- 
tinuum spectra of each source. However, source fluxes at 
100 GHz should be interpreted carefully due to the pos- 
sible contribution of free-free emission at this frequency. 
For QUaD sources with an IRAS counterpart, IRAS far in- 
frared (FIR) fluxes can be used to identify ultracompa ct H II 
(UCH II) re gions using the Wood-Church weU (WC; ' Wood| 
& Churchwell 1989a) criterion: UCH II regions are ionized 



by O stars, which have very similar flux density distributions 
from object-to-object ( [Wood & Churchwell] |1989b[ ) — the 
distribution of UCH II sources in the FIR color-color plane 
should therefore be tightly restricted, an observation which 
is the basis of the WC criteria of \og{F(,o / Fn) > 1.30 and 
log(F25/^)2) > 0.57, where F\ is the wavelength of the IRAS 
band in microns. The QUaD catalog includes a field indi- 
cating whether or not each source satisfies the WC criteria; 
if free-free emission is important, the shape of the flux den- 
sity distribution will be distorted away from that expected 
of a UCH II region and indicate significant free-free emis- 
sion. However, since we cannot assume that all sources in 
the QUaD catalog are UCH II regions, we prefer to use the 
WC criteria simply to indicate whether or not a QUaD source 
is a UCH II region. Of the 505 unique sources in the / catalog, 
182 satisfy the WC criterion, or 36% of all sources detected 
by QUaD are UCH II regions. 

4. 1 . Catalog Field Description 

The fields present in the total intensity source catalog (Ta- 
ble [3]) are designated as follows, with major and minor axes, 
intensities and uncertainties tabulated for each frequency 
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band. If a source was spatially matched between frequency 
bands, the 100 GHz coordinates are quoted. 

1: Source ID 

2: RA: Source right ascension (J2000) in de- 
grees. 

3: dec: Source declination (J2000) in degrees. 
4: CTmnj ' Major axis in arcmin (100 GHz). 
5: cr„„„: Minor axis in arcmin (100 GHz). 
6: /: Source flux in Janskys (100 GHz). 
7: SI: Source flux uncertainty in Janskys 

(100 GHz). 
8: (J„,aj: Major axis in arcmin (150 GHz). 
9: (J,„i„: Minor axis in arcmin (150 GHz). 
10: /: Source flux in Janskys (150 GHz). 
11: 61: Source flux uncertainty in Janskys 

(150 GHz). 
12: aj: Source spectral index. 
13: 6aj: Source spectral index uncertainty. 
14: Alt. name: Alternative name; matched to ei- 
ther IRAS-PSC or PMN catalog. 
15: WC: Indicates whether the source satisfies 
the Wood-Churchwell criteria for ultracom- 
pact H II regions. 

The fields present in the polarized intensity source catalog 
(Table |4|i are designated as follows, with total intensity and 
polarized intensity and angle tabulated for both frequencies. 



1: 


Source ID 


2: 


RA: Source right ascension (J2000) in de- 




grees. 


3: 


dec: Source dec (J2000) in degrees. 


4: 


I: total intensity source flux in Janskys 




(100 GHz). 


5: 


61: total intensity source flux uncertainty in 




Janskys (100 GHz). 


6: 


P: Polarized intensity source flux in Janskys 




(100 GHz). 


7: 


6P: Polarized intensity source flux uncer- 




tainty in Janskys (100 GHz). 


8: 


(p: Source polarization angle in degrees 




(100 GHz). 


9: 


6(t>: Source polarization angle uncertainty in 




degrees (100 GHz). 


10: 


I: total intensity source flux in Janskys 




(150 GHz). 


11: 


61: total intensity source flux uncertainty in 




Janskys (150 GHz). 


12: 


P: Polarized intensity source flux in Janskys 




(150 GHz). 


13: 


6P: Polarized intensity source flux uncer- 




tainty in Janskys (150 GHz). 


14: 


(p: Source polarization angle in degrees 




(150 GHz). 


15: 


6(t>: Source polarization angle uncertainty in 




degrees (150 GHz). 


16: 


aj: Total intensity source spectral index. 



17: 6aj: Total intensity source spectral index 
uncertainty. 

18: ap: Polarized intensity source spectral in- 
dex. 

19: 6ap: Polarized intensity source spectral in- 
dex uncertainty. 

20: Alt. name: Alternative name; matched to ei- 
ther IRAS-PSC or PMN catalog. 

The typical error in each position coordinate are calculated 
from the distribution of position uncertainties taken over all 
sources; we find ctj. =1!^'^ arcmin and cr^- =^g'J arcmin at 
100 and 150 GHz respectively. Uncertainties in angular size 
from fitting each source to an elliptical gaussian function are 
~ 0.2 arcmin. We find typical flux uncertainties of ct/ =!l^ 13 Jy 
and (7/ =11)'^^ Jy at 100 and 150 GHz respectively. 

4.2. Source Distribution with Galactic Latitude 

Figure |4] shows the distribution of discrete sources as a 
function of galactic latitude b, after correction for survey cov- 
erage (a smaller range of b is sampled at lower Decl.). The 
median of the distribution is -0.07° and -0.04° at 100 and 
150 GHz respectively. This negative offset is within one 
beamwidth of = at both frequencies, but supports the re- 
sults of other surveys, such as IS chuUer et al. (2009), who 
found the peak of the distribution to be -0.09 , at higher 
(19.2") angular resolution. No obvious explanation for this 
offset is given, though |Schuller et al.| ( |2009| l suggest the 
slightly positive galactic latitude of the Sun as a possible 
cause, or alternatively the presence of molecular clouds which 
obscure IR sources; the latter is discussed in the context of 
sources near the galactic center in |Hinz et al.| ( |2009) l. 

0.81 ^ . ^ =• ^ . ^ 1 

0.7- 
0.6- 



0.5- 




b (deg) 

Fig. 4. — Probability distribution of detected sources in / as a function of 
galactic latitude b. The solid (dashed) line is for sources detected at 100 (150) 
GHz. At both frequencies, the distribution peaks below h = 0. 



4.3. Source Counts 

Figure |5] shows the differential source counts of the catalog 
as a function of total intensity flux S. Fitting to a power-law 
distribution dN/dS oc S''" in the range 10 < 5 < 300 Jy, we 
find 75 100 = -1 .8 ± 0.4 at 100 GHz, and 75 150 = -2.2 ± 0.4 at 
150 GHz. 

If dust dominates the millimeter/sub-mm source emission, 
their fluxes are proportional to the masses of star-forming 
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Fig. 5. — Differential source counts from the survey as a function of flux. 
Black is for 100 GHz data, gray is 150 GHz. The best-fit slopes are shown 
in the plots as dot-dashed lines with the same color coding; numerical values 
and uncertainties are shown in the lower left of the plot. 

cores M (e.g. Enoch et al. 2006), and the slope of dN/dM can 
be used to constrain the slope of the IMF. Four caveats prevent 
conversion of QUaD source fluxes to core masses. First, the 
resolution required to observe individual cores (as opposed 
to clumps) is approximately 30", a factor ^ 10 higher than 
the QUaD 150 GHz band. QUaD sources could in principle 
contain more than one core, biasing the measurement of core 
masses. Second, since the 100 GHz band flux could contain 
a substantial free-free contribution, calculating masses at this 
frequency is not possible without further information on the 
relative contribution of free-free. This is less of an issue at 
150 GHz, where the QUaD data should be dominated dust. 
Third, the mass conversion also requires a distance estimate 
to the core which are not readily available for each source in 
the catalog. Fourth, the measurement of the slope may be 
subject to systematic error due to the presence of the diffuse 
background. Simulations in Appendix |A.4.2| show that if a 
diffuse background is present the slope of dN/dS may not 
be well-described by a single power-law. However, this effect 
depends on the model used for the diffuse emission in the sim- 
ulations, namely the amplitude of the background relative to 
the sources, and the power in diffuse substructure. Caution is 
thus advised when interpreting the slope results quoted above, 
though Figure[5]indicates that power-law behavior is observed 
above 10 Jy, and thus the contribution of diffuse emission in 
this flux regime is not important. 

Due to these caveats, we caution against overinterpretation 
of the measured slope of dN/dS, and refrain from assigning a 
mass to each source and from converting the slope of dN/dS 
to the slope of the IMF. 

4.4. Spectral Index Distribution 

The spectral index distribution in total intensity, Pr(ai), is 
computed following Muchovej et al.| ( [20 10V For each source 
/', the spectral index probability distribution Prjiai) is calcu- 
lated by generating flux distributions at each frequency from 
the central value and noise distributions, and then combining 
the flux distributions. The spectral index distribution for the 
sample is then the normalized sum of the Prj(aj), i.e. 



Pr(a,) = 



(3) 



Figure |6] shows Priai) for sources matched between the 
QUaD bands. 




Fig. 6. — Source total intensity spectral index distribution Pr{aj). 

The spectral index of the sources are somewhat flatter than 
those found at higher frequencies (e.g. Desert et al. 2008), 
peaking at a/ ~ 0.25; this could be due to the contribution 
from other emission components at 100 GHz, raising the flux 
at this frequency above that expected from dust alone and 
therefore flattening the spectral index. Simulations indicate 
that the spectral index distribution can be slightly skewed 
towards larger aj by background contamination and source 
confusion (see Appendix A. 4. 3 1; the center of the distribu- 
tion shifts by ^ 0.2 in aj. Finally, there is evidence of un- 
accounted emission processes at 100 GHz in the QUaD data 
(see Map Paper), which would also shift the spectral index 
distribution to lower values; detected sources are generally 
not faint enough that flux boosting is important. 

Our analysis does not account for flux boosting due to noise 
on account of the larger (systematic) effect of contamination 
due to the diffuse background and map filtering effects. Rig- 
orous Bayesian methods to determine the spectral index dis- 
tribution of sources exist in the literature (e.g. [Crawford et JI] 
|2010| [Vieira et aL]|2010| ), but do not account for the effect 
of an unknown background, which is the largest contami- 
nant to source fluxes in the galaxy as demonstrated in Ap- 
pendix A.4. 1 We therefore do not pursue such an approach; 



the increased frequency coverage of current-generation satel- 
lite experiments such as Planck and Herschel may allow an 
improved treatment of the diffuse background, enhancing the 
extraction of discrete galactic sources and their spectral in- 
dices. 

4.5. Source Clustering 

Figure [T| demonstrates that source locations in the QUaD 
survey are highly correlated. To quantify source clustering, 
we construct the two-point angular correlation function w{9), 
defined as the excess probability of finding a source within 
angle 9 ± A9 of another source, Hj(9), compared to the same 
probability in a distribution of sources with random spatial 
positions, Hr{6): 

w(6')=-f^-l. (4) 



Hr(9) 



Enoch et al. ( 2006 1 model w as a power law in units of pro- 
jected physical separation r, w(r) oc r'^', and use the slope as a 
method of comparing the spatial properties of cores in differ- 
ent molecular clouds; the authors suggest that different slopes 
may provide insight into the processes dominating core frag- 
mentation. In the QUaD survey, we compute the correlation 
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function in angular units w{9) cx 9'^", by constructing Hci(O) 
from the data, and HriO) from one realization of Siml (a sim- 
ulation with sources distributed randomly over the QUaD sur- 
vey). Differing survey areas at each frequency and variations 
in survey sensitivity are then accounted for 

Any survey over a large range of galactic latitude faces the 
problem that the distribution of sources is anisotropic, with 
w(d) poorly defined at large galactic latitudes due to the shape 
of the galaxy projected on the sky. In addition, at much 
smaller separations w(9) is not well reconstructed due to the 
large probability that a neighbouring source has a low flux, as- 
suming a power law source count dN/dS oc S''\ with 75 < 0. 
Therefore, though a bright source will be detected in the sur- 
vey, its fainter neighbour is likely to lie below the noise or 
confusion limit, preventing accurate reconstruction of w{6) at 
small 9. Simulations indicate that w(9) is well-recovered in 
the range 0.4° < 9 <2°, (see Appendix |A.4.4| i and therefore 
these limits are used to fit a power-law to the correlation func- 
tion of the data. Only sources from the fourth quadrant are 
used, because the gap in survey coverage between the third 
and fourth quadrants introduces artefacts into w; since sources 
in the third quadrant account for < 20% of all sources at each 
frequency, the calculated slope is not affected by removal of 
these sources. 

Figure [7] shows the results. We find a power-law slope 
of 70.100 = -1 .21 ± 0.04 and 761.150 = ~1 -25 ± 0.04, consistent 



with the value found by Enoch et al. 



pOOel ) for the Bolocam 
observations of the Perseus molecular cloud, w(r) oc r"'^^. 
Since a single distance is assumed to Perseus, there is a one- 
one mapping between r and 9, implying that the same cor- 
relation function slope applies to sources on large and small 
angular scales (QUaD and Bolocam respectively). However, 
it should be noted that due to differing resolutions, Bolocam 
and QUaD measure the angular correlation function of differ- 
ent types of source; cores in the case of Bolocam, and clumps 
in the case of QUaD. Therefore, while the correlation func- 
tion slopes are consistent, it is not clear that the correlation 
functions measured by each experiment are directly related, 
and thus caution is advised when comparing these results. 
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Fig. 7. — Angular correlation function w(9) for sources in QUaD survey; 
Blue is for 100 GHz data, red is 150 GHz. Only fourth quadrant data is 
used because the gap in survey coverage between the third and fourth quad- 
rants introduces a discontinuous range of source galacic longitudes. Since the 
number of sources in the third quadrant is only a small fraction of the total 
survey, the effect on the analysis is small. 



4.6. Polarized Sources 



Polarized sources are of particular interest due to their 
scarcity and the fact that they offer a means to probe small- 
scale magnetic field structure in the galaxy. Maps of de- 
tected polarized sources are examined visually in order to re- 
ject beam-scale optical effects (such beam offsets and/or dif- 
fering beam ellipticities between two PSBs within a feed). 
These can cause spurious source detection at fractional po- 
larization (^ 1 % or less) — see Instrument Paper for further 
details. Figure [8] shows the 150 GHz U mnp of RCW 38 and 
a simulated point source of low fractional polarization. The 
'quadrupole' polarization pattern is observed in the two cases, 
implying that the apparent polarization of RCW 38 is an in- 
strumental effect rather than real polarized signal. Polarized 
sources exhibiting such a pattern are visually rejected from 
the catalog. 



150 GHz U RCW 38 



150 GHz U Sim 





1. 







5 134.6 
RA (dcg) 



230 229.6 
RA (dcg) 



Fig. 8. — Illustration of optical effect for low (< 1%) fractional polarization 
sources. Left is 150 GHz U map centred on RCW 38; right is a 150 GHz U 
map of a simulated source showing a similar optical effect. The colorscale on 
each plot is the same. 

Properties of remaining polarized sources in the QUaD sur- 



vey are presented in Table B 
polarized intensity in Figure 
in more detail. 



with images of each in total and 
91 Below we discuss each source 



4.6.1. 284.33-0.36: RCW 49 

RCW 49 is a bright H H region covering 90' x 70', which 
is being ionized by the rich, compact star cluster Westerlund 
2 (e.g. Furukawa et al. 2009). Numerous total intensity ob- 
servations exist from the radio to X-ray wavebands; however, 
existing studies of this sourc e in polarization near t he QUaD 
bands have been restricted to |Dickinson et al.| ( 2007 ), who ob- 
served several southern H n regions with the CBI telescope at 
31 GHz. Their measurements of RCW 49, at 6.78 arcmin 
angular resolution, provide an upper limit on the 31 GHz 
polarization fraction of 0.24%, limited by instrumental leak- 
age from Stokes I to Q and U. Below 31 GHz, the RCW 
49 emission is dominated by free-free, as indicated by the 
total intensity spectral i ndex a/.j?cw49 = ~ 0.220 ± 0.074 be- 
tween 2.7 and 15 GHz d Dickinson elaLjllOOTl l. The QUaD 
total intensity counterpart to RCW 49, 284.33-0.36, also in- 
dicates a flat spectral index between 100 and 150 GHz of 
a/,RCH'49.e(7aD = -0.08 ±-0.01. Polarizcd emission is detected 
in the QUaD data at 100 GHz with a polarization fraction 
of 0.019 ±0.0076 and the polarization vectors aligned pre- 
dominantly east-west (see Figure|9]l. The absence of detected 
polarization at 150 GHz indicates the emission may not be 
thermal in nature; this idea is supported by the total intensity 
spectral index, though this measurement is likely biased flat 
by the presence of free-free emission. While free-free is not 
intrinsically polarized, it may cause polarization by Thom- 
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interferometric instruments at comparable frequencies. This 
result is not surprising, given the fa ctor ~ 20 lower resolu - 

||2006l), 

and a factor ' ' " ' 



tion compared to BIMA at 3.5mm (Macquart et al. 

6 lower than SCUBA at 2mm ( |Aitken et ai.| 



Fig. 9. — Total intensity images withnolarization vectors overlaid on each source detected in polarized intensity — source fluxes, spectral indices and alternative 
identifications are presented in Table|4] Top row is 100 GHz, bottom row is 150 GHz, only polarization vectors with signal-to-noise > 5 are plotted. From left 
to right: 284.33-0.36 (RCW 49), 359^93-0.06 (Sagittarius A*), 0.18-0.06 (Galactic Center Arc), and 0.63-0.05 (Sagittaiius B2). For 284.33-0.36, the contours 
run from 1.6 to 16.1 MJy/sr in steps of 1.6 MJy/sr at 100 GHz, and from 2.3 to 23.3 MJy/sr in steps of 2.3 MJy/sr at 150 GHz. For the remaining sources, the 
contours run from 2 to 20 MJy/sr in steps of 2 MJy/sr at 100 GHz, and from 5 to 50 MJy/sr in steps of 5 MJy/sr at 150 GHz. 

son scattering at the edges of the H II region, resulting in 
tangentially polarized radiation at the cloud edges. Since the 
QUaD 100 GHz polarization vectors are largely aligned over 
the source area, free-free polarization at the cloud edges can 
be ruled out. Synchrotron radiation is a further possibility, but 
unlikely given the physical nature of the source. It is therefore 
possible that instrumental effects other than those illustrated 
in Figure[8]are present. 

4.6.2. 359.93-0.06: Sagittarius A* 

High frequency polarized observations of Sagittarius A* 
(Sgr A*, detected in the QUaD survey as 359.93-0.06) are 
important due to the constraints they provide on relativistic 
jets and accretion processes in black holes (e.g. Agol||200"0 
IQuataert & Gruzinov||2000l |Melia et al.||20oT] i. The rota- 
tion measure (RM), which is proportional to the electron den- 
sity and magnetic field component integrated along the line 
of sight, can be determined from multi-frequency observa- 
tions of the linear polarization fraction — measurements of 
the RM constrain the mass accretion rate of the black hole and 
thus rule out certain classes of accretion model. Observations 
of linear po larization in Sgr A*, first detected with SCUBA 
by Aitken et al. ( 2000 ), have been studied with interferometric 
imaging at frequencies above ~ 100 GH z, typical ly with res- 
olutio n < 10" or be tter (Macquart et al. 2006; M arrone et al. 
|2007| l. At 3.5mm, M acquart et al.| ( |2006) ) find a fractional lin 



2000). At 3.5' angular resolution at 150 GHz, the polai'ized 
emission from Sgr A* is smeared out and the polarization 
fraction therefore reduced due to the contribution from the 
diffuse background or unpolarized sources within a QUaD 
150 GHz beam of Sgr A*. Though QUaD observations span 
four months and could constrain the variability of the polar- 
ized emission (a useful diagnostic of processes intrinsic to the 
source), at such low resolution and signal-to-noise it is not 
possible to measure this effect. The apparent low polarization 
fraction of this source may result in false detection due to an 
optical effect in the telescope, similarly to 284.33-0.36, and 
should therefore be treated with caution. 

4.6.3. 0.18-0.06: The Galactic Center Arc 

0.18-0.06 is clearly associated with GO.2-0.0, also known 
as the Galactic Center Arc. This source has an extent of 
^ 25 arcmin along its long axis perpendicular to the galactic 
plane, is approximately symmetric with respect to the galactic 
equator, and is among the brighter sources close to the galac - 
tic center (e.g.'Alten hoff et al.|1979HLis & Carls troml 1994 



ear polarization of2.1%±0.1%, while larger values are found 
at higher frequ encies (~ 5% and ~ 9% at 230 and 340 GHz 
respectively, Marrone et"aL] ( |2007| )). [Aitken et aT] ( |2000 1 
measured an observed polarization fraction of 2.9% ± 0.3% 
at 2mm with 33.5" resolution. 

In the QUaD data, Sgr A* is detected in polarization at 
150 GHz only (see Figure [9] and Table [4]), but the polarized 
flux is unconstrained; it appears that the polarization fraction 
is < 1% at 150 GHz, considerably less than that found by 



High resolution observations (e.g. |Yusef-Zadeh et al.||2004| l 
indicate this filamentary structure contains the largest concen- 
tration of nonthermal radio filaments in the galaxy. Though 
this object appears extended and unassociated with any dis- 
crete source in the QUaD / data (see Figure |9] and Table |4]), 
the radio filaments cross three bright H II regions unresolved 
by QUaD (e.g. |Yusef-Zadeh et al.|[T986l [Reich et al.|[2000l ): 
GO.16-0.15, GO.18-0.04 and G0.1h-0.08. Using the Green 



Bank Telescope, Law et al. (2008) determined that this radio 
arc has a non-thermal spectrum of -0.54 ± 0.09 between 4.85 
and 8.5 GHz (resolution 2.5 and 1 arcminute respectively), in 
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support of the idea that the emission is not from cold dust, but 
rather from either monoenergetic el ectrons or an elec tron dis- 
tribution with a low energy cutoff (Reich et al. 120001). Obser- 
vations at 4.8 GHz using the VLA (Yusef-Zadeh et al. 1986 ) 
demonstrate that while thermal emission dominates GO. 18- 
0.04 and G0.1h-0.08, the low-frequency nonthermal polarized 
emission in the arc is primarly due to the H II region GO. 16- 
0.15; the high 30%) fractional polarization observed pro- 
vides further evidence that the polarized emission is due to 
synchrotron. 

We find a peak polarization fraction of ^ 10% at both 100 
and 150 GHz, and a polarized spectral index of -1.04 ±0.17, 
indicative of a source dominated by synchrotron emission. 
Since we do not detect a discrete total intensity source as- 
sociated with the polarized arc, we estimate the polarization 
fraction from the raw / images and the polarized flux mea- 
sured by the source extraction algorithm. Due to its spatial 
extent above the beam scale in the QUaD survey, it is doubt- 
ful that the radio arc's polarized emission at these frequencies 
can be solely attributed to GO. 16-0.15, which would be likely 
be unresolved. The high level of uniformity of the polariza- 
tion seen in Figure|9]indicates a highly ordered magnetic field 
which must exist over the full extent of the arc, rather than 
localized to a single H II source. The polarized vectors are 
largely aligned parallel to the plane of the galaxy, almost per- 
pendicular to the polarization from diffuse emission (see Map 
Paper), indicating a strong local deviation from the galactic 
magnetic field. 

4.6.4. 0.63-0.05: Sagittarius B2 

The polarization of Sagittarius B2 (Sgr B2) has been well- 
studied in the radio, sub-millimeter and far-infrared (e.g. 
Greaves et al.|[T995l |DoweU||1997[jNovak et al.||1997[ [Eo^ 



ell et al.||1998| [Jones et al.||2010| l; the QUaD data fill in 
the millimeter portion of the spectrum of this giant molecu- 
lar cloud. Observations of linearly polarized emission from 
magnetically aligned dust grains in such clouds can be used 
to determine the orientation of the local magnetic field. In 
the sub-millimeter, where the dust is optically thin, the po- 
larization is due to emission of the dust grains preferentially 
along their long axis; [Greav es et al.| ( [T995) observed Sg B2 
at 800/im and 30" resolution, finding polarization fractions 
in the range 0.8% -2.6% with the polarization orientation 
approximately north-south. Polarization observations at far- 
infrared wavelengths include both emission and absorption 
effects. The former is due to dust grain emission similarly 
probed by sub-millimeter observations and results in polariza- 
tion vectors aligned perpendicular to the local magnetic field 
direction, while the latter is caused by absorption from cold, 
magnetically aligned dust grains in regions of high optical 
depth, with correspondi ng polarization vect ors aligned paral- 
lel to the magnetic field. Novak et al. ( 1997| demonstrate that 
at 1 15 fim with a 35" beam, Sgr B2 is resolved into 'core' and 
'envelope' regions, with polarization in the core and envelope 
dominated by absorption and emission respectively. In their 
study, the envelope fractional polarization ranges from 2-A%. 
At 350 /im, the effects of absorption are diminished due to 
the increasing co ntribution of dust em ission and the decreas- 
ing optical depth; Dowell et al. ( 1998 1 show this observation- 
ally, with the core polarization (^ 1 % fractional polarization) 
smaller than that in the envelope (^ 2.8%) at 20" resolution. 

Based on the above considerations, in the QUaD bands the 
polarization of Sgr B2 should be dominated by dust grain 
emission processes. The best measurement of polarization 



fraction is at 150 GHz, where we find a value of 1% ±0.2%. 
It is not surprising that the 150 GHz polarization fraction is 
smaller than at 350 /im. Since the QUaD beam at this fre- 
quency is a factor 3.5 x 60/20 ~ 1 1 larger, the core and enve- 
lope are not resolved into separate components, and thus the 
higher envelope polarization fraction is biased low by the less 
strongly polarized core region. The east-west orientation of 
the 150 GHz QUaD polarization vectors (see Figure |9]), sim- 
ilar to those at 350 /im (Dowell et al. 1998 ), further indicates 
that we are observing polarized dust emission from the enve- 
lope. 

4.6.5. An Extended Polarized Source: Further Filtering 
Considerations 

For the sake of simplicity, the same aggressive background 
filtering (small dhck relative to (Jbeam) is used in both total and 
polarized intensity. This choice was motivated by the bright 
diffuse emission in total intensity: a larger ahck would result 
in a smaller systematic loss of signal from point sources but 
an increased contribution from the background — see Ap- 
pendix |A.5| However, this choice may be relaxed in polar- 
ized intensity because the polarized diffuse emission is rela- 
tively faint (polarization fraction < 2%, see Map Maper). The 
present catalog is therefore not optimized for the detection of 
polarized sources, and more sources could in principle be ex- 
tracted from the maps. 



As an illustration. Figure 10 shows total and polarized in- 
tensity maps for a visually identified cloud with appreciable 
polarization at 150 GHz (up to 10% fractional polarization). 
Since this source has a relatively large extent in comparision 
to the beam scale, it is undetected by the source extraction al- 
gorithm because most of the polarized flux is filtered out in 
the background removal step. 

This object is host to several discrete sources, namely 
344.99h-1.79, 345.01+1.53, 345.37+1.42 (resolved into 
345.38+1.41 and 345.49+1.46 at 150 GHz), 345.22+1.03, 
and 344.96+1.23 (detected at 150 GHz only). Sources 
345.01+1.53 and 345.08+1.59, detected at 100 and 150 GHz 
respectively, are cleary associated with the same source, but 
have observed centroids separated by more than the 0.02° 
required to meet the spatial matching criterion. The polar- 
ized emission in this region is only partially correlated with 
the positions of the discrete / sources (see Figure [T0| , and 
the polarization vectors are largely oriented perpendicular to 
those from the bulk polarized galactic emission. This indi- 
cates that the diffuse dust is subject to a local magnetic field 
strong enough to overcome that of the galaxy as a whole, de- 
spite its small galactic latitude of ^ 1.4°. Clearly this source 
is worthy of detailed follow-up study. 

5. DISCUSSION AND CONCLUSIONS 

We present a catalog of discrete sources extracted from the 
QUaD galactic plane survey, which spans approximately 245- 
295° and 315-5° in galactic longitude I and -4 to +4° in galac- 
tic latitude b — a total of ^ 8()0 square degrees coverage in 
Stokes /, Q, and U at 100 and 150 GHz, with resolution 5 
and 3.5 arcminutes respectively. Simulations of a toy model 
galaxy including spatially clustered point sources and diffuse 
emission indicate a 90% completeness flux of 5.9 (2.9) Jy at 
100 (150) GHz in /, and 20.3 (1.1) Jy in polarization at 100 
(150) GHz. The high 100 GHz completeness in polarization 
is due to source confusion in the larger beam, and is depen- 
dent on the parameters used for the toy galactic model; ran- 
domly distributed sources yield a completeness of 1.3 Jy — 
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1 344.99+1.79 

2 345.01+1.53 




3 344.96+1.23 

4 345.37+1.42 

5 345.49+1 

6 345.22+1.03 5" 




Fig. 10.— Large-scale polarized cloud at 100 GHz (left) and 150 GHz 
(right) undetected by source extraction algorithm — only polarization vectors 
with signal-to-noise > 4 are plotted (black lines). At 100 GHz, the contour 
scale (gray lines) runs from 0.32 to 3.2 MJy/sr in steps of 0.32 MJy/sr, and 
from 0.77 to 7.7 MJy/sr in steps of 0.77 MJy/sr at 150 GHz. Black squares 
indicate the locations of discrete sources detected in the / maps, with the 
names of the sources as indicated at upper left of each plot. 



see Appendix A. 2 for a full discussion. At a signal-to-noise 



threshold of 5 (3) in total (polarized) intensity, the catalog is 
98% pure in / at both frequencies, and 97% (92%) pure in 
polarization at 100 (150) GHz. Simulations without a dif- 
fuse background are used in the total intensity computation 
because substructure in diffuse emission, detected as discrete 
sources, biases the purity low. While this could also affect 
the real catalog, the high percentage of IRAS-PSC counter- 
parts to QUaD sources (97% and 87% at 100 and 150 GHz 
respectively) indicates that this effect is likely very small. 
The polarized diffuse background, with fractional polariza- 
tion ~ 2% (see Map Paper), is faint enough that it does not 
bias the purity of the catalog at the signal-to-noise threshold 
of the survey. Instrumental effects prevent detection of polar- 
ized sources with polarization fraction ^ 1% or less. 

In total intensity the catalog contains 505 unique sources, 
of which 239 are spatially matched between frequency bands, 
with 50 (216) detected at 100 (150) GHz alone. 

The / flux distributions are well-approximated by a power 
law over more than two orders of magnitude above ~ 10 Jy 
at both frequencies. We find power-law slopes of 75 100 = 
-1.8 ±0.4 at 100 GHz, and 75,150 = -2.2 ± 0.4 at 150 GHz; 



the latter is consistent with Rosolowsky et al. ( 2010), who find 
-2.4 ± 0. 1 at 268 GHz with Bolocam at higher resolution. The 
flatter slope at 100 GHz may be the result of resolution effects 
due to the larger beam at this frequency. Simulations indicate 
that if the diffuse background contributes spurious sources, as 
expected the recovered source flux distribution does not ac- 
curately follow the underlying distribution; however, as dis- 
cussed above, the high percentage of QUaD sources spatially 
matched to IRAS indicates this effect is insignificant. 

The spectral index probability distribution of sources in to- 
tal intensity is found to peak at a 0.25, flatter than ex- 
pected for sources whose emission is dominated by ther- 
mal dust. Simulations indicate the diffuse background does 
not strongly influence source spectral indices; the flatness is 
therefore likely due to free-free emission, which becomes sig- 
nificant at ^100 GHz and below. At this frequency, free-free 
results in higher fluxes than expected from dust alone, shifting 
the spectral index distribution to lower values. 

We explore the clustering of galactic sources by fitting the 
two-point correlation function to a power-law using the / 
source locations. Simulations indicate that the underlying cor- 
relation function slope can be accurately reconstructed in the 



range 0.4° < 6 <2°, with 9 the angular separation between a 
pair of sources. The correlation function breaks down at larger 
angular scales because so few (< 1%) of sources are located 
beyond \b\ = 3°. At galactic latitudes smaller than 0.4°, w(d) 
is not well reconstructed because for a power-law dN/dS oc 
S] with 75 < 0, the survey does not detect most neighbors of 
a source bright enough to be included — one must extend the 
search to large angular separations before enough neighbors 
are detected for accurate reconstruction. Fitting to the QUaD 
/ catalog data in the range 0.4° < 9 < 2° , we find power-law 
slopes of 7e,ioo = -1.21 ±0.04 and 79,150 = -1.25 ±0.04 at 
100 and 150 G Hz respectively. The se are consistent with the 
value found by 'Enoch et al. ( 2006 1 Bolocam observations of 



the Perseus molecular cloud, w(9) oc 9~^ '^^, though the results 
are not directly comparable on account of the different sources 
probed by QUaD (clumps) and Bolocam (cores) due to their 
differing angular resolution. 

97% (87%) of the sources detected at 100 (150) GHz have 
IRAS-PSC counterparts. These large fractions indicate that 
most of the clumps detected in the survey are past the prestel- 
lar phase and have envelopes heated by protostars. This ob- 
servation might be expected, given that the QUaD frequency 
bands lie far from the spectral peak: Only these sources are 
bright enough in the Rayleigh-Jeans portion of the spectrum 
to be detected in the QUaD survey, unlike prestellar or star- 
less sources. Since the QUaD survey is sensitive to free- 
free emission as well as dust, particulary in the 100 GHz 
band, sources might also be detected if their gas is sufficiently 
ionized to produce free-free but their envelopes are yet to 
thermalize. However, the small fraction of 100 GHz QUaD 
sources unmatched to IRAS-PSC imply that almost all the 
detected sources at this frequency do have a thermal compo- 
nent. At 150 GHz, the larger unmatched fraction is likely due 
to single IRAS sources being resolved into two sub-clumps 
by the the higher QUaD resolution at this frequency — only 
one of these sub-clumps can be spatially matched to the IRAS 
source. 

Of the sources with an IRAS counterpart, 182 satisfy 
the Wood-Churchwell criteria for ultracompact H II re- 
gions (Woo d & Churchwell||1989"a] ), providing new spectral 
constraints on this class of object. 

Four compact polarized sources were detected by the auto- 
mated source-finding algorithm: 284.33-0.36 (IRAS 10227- 
5730 or RCW 49), 359.93-0.06 (Sagittarius A*), 0.18-0.06 
(IRAS 17431-2846 or Galactic Center Arc), and 0.63-0.05 
(IRAS 17440-2823 or Sagittarius B2). One additional ex- 
tended source was located 'by eye' from the raw Q and 
U maps; this object appears host to several discrete to- 
tal intensity sources, including 344.99h-1.79, 345.01h-1.53, 
345.37+1.42, and 345.22+1.03. The brightest polarized 
source is 0.18-0.06, which does not have an obvious discrete 
counterpart in total intensity, but has a polarization fraction 
of ^ 10% if the diffuse background is used as a measure 
of /. It has a polarized flux of 7.91 ±0.33 (4.90 ±0.32) Jy 
at 100 (150) GHz, and a polarized spectral index of ap = 
-1.04 ±0.17, indicating a synchrotron emission source. Its 
detection against a polarized background implies that there is 
a strong local deviation from the galactic magnetic field. 

Less than 1 % of the sources detected in / have a polarized 
counterpart. If discrete sources do not harbor strong local 
magnetic fields or shielding, dust grains in their envelopes 
will align with the large-scale galactic field. The only way 
to separate diffuse from discrete polarized emission would 
then be via morphology (similar to /) or spectrally, since the 
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orientation of the polarization would be similar for diffuse 
and discrete sources. Alternatively, the discrete total inten- 
sity sources may have fractional polarization < 1%, as might 
be expected from a star-forming clump, in which case instru- 
mental effects prevented detection of their polarized emission 
here. Discrete sources may therefore not be a significant con- 
tributor to the low-latitude galactic polarized emission. More 
sensitive observations (such as from the Planck satellite) will 
be needed to better study the polarization of these sources, 
and the role of magnetic fields in star-forming regions. 

The QUaD catalog may prove useful for a variety of ad- 
ditional purposes. Total intensity source fluxes could better 
measure the continuum spectra of clumps in conjunction with 
independent data sets, improving the separation of different 
emission components and tightening constraints on dust emis- 
sivity and gas temperatures. The maps provide upper limits to 
source polarization, allowing a statistical study of polarized 
contribution to anomalous emission similar to lDickinson et al.l 
([2007 ). Finally, the catalog provides a cross-check of astrom- 
etry and absolute calibration for instruments with access to 
the southern hemisphere. 
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APPENDIX 
A. SIMULATIONS 

The source extraction algorithm presented in Section|3]is tested using simulated distributions of point sources, with and without 
a toy model diffuse background. These simulations are used to determine the survey completeness, purity, the accuracy of 
recovered individual source parameters, the accuracy of recovered source distribution parameters, and the effect of choice of 
background smoothing kernel width ahck- 



A.l. Galactic Model Generation 



The methods described in Appendices A. 1.1 and A. 1.2 below are used to generate four types of galaxy simulation: random 
and correlated spatial distributions of point sources, with and without a diffuse background component — Table |2] summarizes 
the properties of each simulation. 



A. 1 . 1 . Simulated Source Populations 

Point sources are placed over the area of pixels occupied by the QUaD survey, using two methods to generate their spatial 
distribution. The first is a simple random distribution, denoted Siml, with a partner simulation (Sim2) also constructed from the 
same spatial distribution of sources but with the inclusion of a diffuse background component (see Appendix |A.1.2| below). The 
second more closely matches the observed clustering of discrete sources: A two-point angular correlation function w (6) = k9~^^ is 
used to model the clustering of sources, and a power-law probability distribution function in galactic latitude b is simultaneously 
implemented to capture the observation that sources tend to be concentrated towards the galactic plane, i.e. p(b)Ab oc b~^db. This 
simulation is called Sim3, and a further simulation, Sim4, is generated by taking the same point source population and adding 
a diffuse component. The power-law exponents and /3 are chosen such that the spatial distribution of sources qualitatively 
matches that in the QUaD data. Since it contains both a diffuse background and spatially clustered point sources, Sim4 is the 
model which most closely resembles the real data. 

The physical properties of the sources are defined as follows: 100 GHz total intensity source fluxes are drawn from a power-law 
model for the source counts (dN/dS oc with 75 = -1.5) between 0.1 and 250 Jy; the normalization is chosen to match the 
average source density (i.e. number per square degree) in the QUaD / catalog. Spectral indices between 100 and 150 GHz are 
generated using a gaussian probability distribution function (p.d.f) of zero mean and unit width; 150 GHz fluxes are generated by 
combining the spectral index and 100 GHz flux. The polarization fraction for each source is a random number drawn uniformly 
between and 20%, while the polarization angle is also a uniform random number between and 180°. 



A. 1 .2. Simulated Dijfuse Background 

Adding a model diffuse background allows its effect on recovered source properties to be assessed. This component is modelled 
using a weighted sum of a single point source map smoothed to different resolutions. 



Point source locations are generated using the correlation function approach described in Appendix A. 1 . 1 for Sim3 and Sim4, 
with a source density ^ 20 times higher than the real data. All sources have the same 100 GHz fluxes Sbck, a spectral index of 
2.18 (the mean spectral index of the diffuse emission as calculated in the Map Paper), and a polarization fraction of 2% in pure 
+Q in galactic coordinates; these sources are placed within the boundaries of the survey in /, Q and U maps at each frequency v, 
giving a set of template maps m,^_,-, where / refers to each Stokes parameter. 

To construct the model of diffuse emission, each m,„ i is first smoothed with a set of circular gaussian kernels of differing 
width aj to yield smoothed maps m,j,,ij, with -1.4 < logi()(f7j) < 0.2 in intervals of Slogioiaj) = 0.2. The smoothed maps are 
then coadded with a different weight wj for each smoothing kernel, with wj oc f^'^, yielding the diffuse model map inhck,u.i' 

mhck.v.i = ^Wjin,^^jj. (Al) 

J 

The Wj amplitudes and Sbck were chosen by requiring that the flux integrated over the whole simulated map matched that in the 
data within ^ 10%. A further consideration is that the relative weights between the wj, and the exponent (defined somewhat 
arbitrarily here to be 7/4), should be chosen such that substructure on different angular scales in the model matches the data; a 
power spectrum analysis may represent the best way of determining the wj but is beyond the scope of this paper. Here, we merely 
note that with the choice of (Tj and wj used above, the diffuse component of the simulations bear a qualitative resemblance to the 



data, as may be seen in Figure 1 1 



A. 1.3. Simulated Maps 

Having determined the positions and physical properties of the sources, they are then placed in a 'source map', mj„-, of 
the same size and pixel resolution as the QUaD maps. A similar map nibck is generated for the diffuse background following 
Appendix |A.1.2[ nisrc alone is the input sky for Siml and Sim3, with nisrc + mhck used for Sim2 and Sim4. For each type 
of point source spatial distribution (random or correlated), the coordinates of each source are the same with and without a 
diffuse component present. This allows the effect of the background on source fluxes to be investigated separately from spatially 
correlated source positions. Note that a different realization of the diffuse component is used for the random and correlated source 
simulations, as may be seen in Figure [TT] 

Simulated detector timestream is interpolated from these maps using the pointing information from each day of real QUaD 
data, and realistic noise added as described in the Map Paper. The simulated signal-nnoise data is then subjected to the same 
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Fig. 1 1. — Comparison of QUaD 100 GHz map to one realization of the simulations described in Appendix |A.1.3| with all maps transformed to galactic 
coordinates. The color scale is MJy/sr. From top to bottom: QUaD data, random point sources only, random sources plus diffuse background, spatially correlated 
sources, and spatially correlated sources plus diffuse background. 



TABLE 2 
Simulation Summary 



Simulation 


Source spatial 


Diffuse 




90% Completeness 






Purity" 






distribution 


emission 


^100 (Jy) 


/i5o (Jy) 


Pm (Jy) 


Piso (Jy) 


P/,100 


PlASO 


Pp. 100 


Pp. 150 


Siml 


random 


no 


1.6 


1.5 


1.2 


0.9 


0.88 


0.89 


0.90 


0.91 


Sim2 


random 


yes 


1.5 


1.5 


1.3 


0.9 


0.43 


0.40 


0.93 


0.91 


Sim3 


correlated 


no 


5.1 


2.4 


13.7 


1.2 


0.98 


0.98 


0.97 


0.92 


Sim4 


correlated 


yes 


5.9 


2.9 


20.3 


1.1 


0.46 


0.43 


0.93 


0.94 



^ Purity estimated at point source signal-to-noise threshold of 5 (3) in total (polarized) intensity. 

filtering and mapmaking steps as in the QUaD data pipeline described in the Map Paper Figure [TT] shows the QUaD 100 GHz / 
map, and maps from one realization of each simulation type. 

The maps are passed through the source extraction algorithm described in Section[3]to generate catalogs of sources in /, Q and 
U at the two QUaD frequencies. To determine quantities such as spectral index and polarization fraction (that is, quantities which 
require the catalog from more than one map), sources are matched using the simple spatial criterion described in Section[3] Many 
sky realizations are processed to build up sufficient statistics to characterize the survey depth and systematic effects. 

A.2. Completeness 



The completeness C(> S), the fraction of input sources recovered above flux S, is shown in Figure 12 for the QUaD frequency 
bands for each set of simulations. The 90% completeness limits (the flux Sgo at which the completeness reaches 90%) for the 
simulated catalogs are summarized in Table [2] 

The 90% completeness limit in total intensity is generally higher than in polarized intensity; this is in part due to the higher 
signal-to-noise detection threshold than that used in polarization (5 compared to 3). The limit in / is also increased due to 
the shape of the completeness curve, where C(> S) increases rapidly to 0.5 in / but then approaches unity more slowly — see 
Figure 12 This effect is attributed to residual 1// noise in the / map, rather than the diffuse background or source confusion, 
since it is present for all simulations. In Siml and Sim2, the 90% completeness limits for / are 1.6 and 1.5 Jy respectively at 
100 GHz, and 1.5 Jy for both simulation types at 150 GHz. The 90% Umit is higher for Sim3, reaching 5.1 (2.4) Jy at 100 (150) 
GHz. Since this simulation contains no background, the effect is solely due to increased source confusion resulting from the 
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Fig. 12. — Completeness estimated from simulated point source distibutions in total (top) and polarized (bottom) intensity. The left pair of columns is for 
100 GHz, right pair for 150 GHz. Within each pair, the left column is for randomly distributed sources with and without a diffuse component, and likewise for 
spatially correlated sources in the right column. Solid black curves correspond to sources with no diffuse background present, dashed black curves are for the 
same source distributions with diffuse emission added. The thin dashed black line indicates 90% completeness. 



spatial clustering of the sources. For Sim4, these numbers increase marginally to 5.9 and 2.9 at 100 and 150 GHz respectively, 
indicating that the spatial clustering of sources is more important than the presence of a diffuse background. 



In polarization, the diffuse background is faint enough that it does not strongly affect the completeness as shown in Figure 12 



A more important effect is the spatial clustering of sources — at 100 GHz, the survey is only 90% complete in polarization by 
13.7 Jy for Sim3 and 20.3 Jy for Sim4. At 150 GHz, these figures drop to 1.2 and 1.1 Jy respectively. This is likely due to three 
effects. First, the survey coverage is smaller at 100 GHz due to the smaller area of the QUaD focal plane at this frequency (the 
100 GHz survey is approximately 7% smaller in polarization than 150 GHz). As a result, some input sources may not lie in 
the 100 GHz survey area; however, for spatially correlated sources, which are clustered close to the plane of the galaxy, this is 
likely a small effect. Second, the larger beam at 100 GHz causes source confusion. The randomly oriented polarization angles 
could result in polarized flux dilution, reducing the completeness. Figure 12 shows that this effect is reduced at 150 GHz, in 
support of this idea. Third, the completeness is a function of polarization fraction and total flux. No source with zero polarization 
fraction, or no faint source with low polarization fraction, will therefore be detected. Since the input source polarization fractions 
lie between and 20%, this effect could also contributed, though it would be expected at both frequencies. In summary, in all 
simulations except 100 GHz Sims 3 and 4, the 90% point source completeness limit in polarization is 0.9-1.3 Jy. 



A.3. Purity 

The 'purity' of the survey p is the number of recovered sources which were matched spatially to the input catalog divided 
by the total number of recovered sources. This is quantified by comparing the input and recovered source catalogs in Sims 1-4. 
Figure 13 shows the purity as a function of signal-to-noise threshold in both total and polarized intensity. The values of p at the 
chosen extraction thresholds (S/N > 5 in /, S/N > 3 in P) are summarized in Table|2] 

Between detection thresholds of 2 < S/N < 8, the purity increases more rapidly with S/N in both total and polarized intensity 
if the spatial distribution of sources is correlated rather than random. This is due to flux boosting of a source by fainter, spatially 
coincident sources which are not resolved themselves. In total intensity, the purity appears poor (40-50%) when the diffuse 
background is present (Sim2 and Sim4), even at a detection threshold of S/N > 12. However, these 'false detections' are not noise 
fluctuations: Investigation of the output catalogs from Sim2 and Sim4 demonstrated that the detected sources unmatched to the 



input catalog were beam-scale or extended sources associated with substructure in the simulated diffuse background. Figure 14 



illustrates this effect for a Sim4 realization, demonstrating how substructure can be erroneously detected as real discrete sources. 
This effect makes 'purity' an ambiguous concept in the context of separating sources from a diffuse background with enough 
power on beam-sized angular scales: If the signal-to-noise of the diffuse background is comparable to that of the discrete sources, 
the purity is not dominated by noise fluctuations, but the inability to distinguish the sources of interest from resolved substructure 
in the diffuse emission. Simulations Sim2 and Sim4 fall into this category, which is why their purity is low even at high signal-to- 
noise thresholds; this effect is simply a consequence of the model parameters chosen in Appendix A. 1.2 As shown in Section|4j 
97% (87%) of the QUaD sources detected in total intensity at 100 (150) GHz have IRAS-PSC counterparts; this indicates that 
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Fig. 13. — Purity as estimated from simulated point source distributions as a function of signal-to-noise (S/N). Tlie left (right) pair of columns correspond to 
purity at 100 (150) GHz. Within each pair, the left plot is for randomly distributed sources with and without a diffuse background (Siml and Sim2), and the 
right plot similar but for spatially correlated sources (Sim3 and Sim4). Survey purity in total (polarized) intensity are plotted in black (gray). The vertical black 
dot-dashed lines coirespond to the S/N ratio used to extract sources in /; the vertical gray dot-dashed line is the same for polarization. 

the QUaD catalog is likely very pure, and that more beam-scale substructure is present in Sim2 and Sim4 than in the real data. 
A power-spectrum approach, advocated in Appendix ] A. 1.2| represents the best way to determine the amount power on different 
angular scales. 




335 340 345 350 



1 (deg) 

Fig. 14. — Illustration of purity in the presence of diffuse background. The image is a subsection of the Sim4 100 GHz / map shown in the bottom panel of 
Figure[l3] Crosses correspond to the locations of input sources, squares are the sources detected from Sim3, and circles are the sources detected in Sim4 (recall 
that the input point source populations are identical for Sim3 and Sim4). In Sim4, false sources are detected from substructure in the diffuse emission — these 
sources have a circle that does not enclose a cross. This systematic effect is the cause of the low purity for simulations which have a diffuse background present, 
as described in Appendix |A.3| 



The diffuse polarized background is faint enough that there is little difference between the purity in each type of simulation. 



Figure 13 shows that a polarized source detection threshold of S/N > 3 results in a catalog that is ^ 90% pure for randomly 



distributed sources, or ^ 100% pure for spatially correlated sources. 

A.4. Source Recovery 

Recovered catalogs from the four sets of simulations are used to estimate how accurately the input parameters of individual 
sources and source distributions can be recovered. While the output distributions for Siml and Sim3 (those without a diffuse 
background) should be insensitive to the choice of input distribution parameters, such as the source counts slope js, the systematic 
biases introduced by the background (simulations Sim2 and Sim4) do depend on the background model parameters 7^, (3, w,, (t,, 
and the overall amplitude. We therefore caution that while the toy model of the diffuse component allows a qualitative impression 
of how source properties can be corrupted, the amount of corruption depends on these parameters to an extent that may differ 
from the real data. 



A.4. 1 . Recovery of Individual Source Parameters 

Figure[T5]compares the recovered source properties to their input values. In total flux S, we show the 16, 50 and 84% percentiles 
of the ratio Sout/Sin, while the same percentiles are shown for polarization fraction difference font- fin and polarization angle 
difference (t>out-<t>in- Both 100 GHz and 150 GHz simulations are shown, though similar behaviour is observed at both frequencies. 

For total intensity S, in all simulations the median S„,i,/Si„ falls as sources get fainter, but is within a few percent of unity down 
to 10 Jy (comparable to the absolute calibration uncertainty of 3.5%). At this flux, the 16 and 84% percentiles of Sout/Sm are 
~ 5% from the median for simulations without a diffuse background (Siml and Sim3). Below 10 Jy, rather than more numerous 
faint sources being boosted to higher fluxes, the systematic defecit in recovered flux is due to filtering of the timestream before 
map coaddition (see Map Paper); this demonstrates that filtering effects are more important than flux boosting due to instrumental 
noise at low fluxes. 

Including diffuse emission results in a wider and more asymmetric distribution of recovered flux due to background contami- 
nation; at 10 Jy in Sim3 (no diffuse emission), the systematic uncertainties are ct- = 0.04 and a+ = 0.07, while the same quantities 
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Fig . 15. — Recovery of simulated point source properties as a function of the input. Left to right, the columns are total intensity source flux S, polarization 
fraction /, and polarization angle (p. S is expressed as a ratio between the output and input values, while / and <p are expressed as the difference of the two 
quantities. Top two rows are 100 GHz, while lower two rows are 150 GHz. Within each pair of rows at a fixed frequency, the simulations used are Siml and Sim2 
in the first row, and Sim3 and Sim4 in the second row. The contours indicate the 16, 50 and 84th percentiles of each quantity as calculated from simulations. Solid 
black contours are for simulations without a diffuse background (Siml and Sim3), while gray dashed contours indicate the simulations with diffuse emission 
(Sim2 and Sim4). 

in Sim4 (including diffuse emission) are ct_ = 0.07 and a+ = 0.36, where ct_ is the difference between the 50th and 16th per- 
centiles of the ratio Som/Sm at 10 Jy, and a+ is the di fferenc e between the 84th and 50th percentiles. The asymmetric errors due 
to background contamination are discussed further in A. 4.2 in the context of the source counts, dN/dS. 

The recovery of point source polarization fraction is not strongly biased in any of the simulations. The median difference 
between input and output polarization fractions is < 1 %, at which point beam systematic effects become important (see Instrument 
Paper for details). The scatter on /„„, generally increases with higher polarization fraction, rising from ~l%at/~l%to 
5% for / = 20%. This is because only the brightest sources tend to be detected at low /, while at high / both bright and faint 
sources are included, increasing the variance of /„,„-/,„. The distribution of /o,„-/,„ is skewed towards positive values; this 
effect is attributed to the addition of noise to the total polarized flux as ap = a-Q + a^. No systematic bias is introduced when 
a diffuse component is present, indicating that the background removal strategy is effective for determining fluxes of polarized 
sources. 

One might ask why the recovered / is not systematically lower than the input, since the polarization angles assigned to simu- 
lated sources are random and therefore should average to zero when sources are confused, as in Sim3 and Sim4? The reason is 
that the simulations use a power-law distribution of fluxes, resulting in many more faint sources per unit solid angle than bright 
sources. In a given resolution element, faint sources will largely cancel each others' polarized flux, while a statistically unlikely 
(but far brighter) polarized source will dominate the polarized emission. 

The simulations show that polarization angle difference 0o,„ -</),„ suffers a systematic shift of < 1°, with 16 and 84 percentiles 
less < 5°. The scatter in (p is not strongly affected by a diffuse component, but increases marginally when sources are spatially 
correlated. This increased scatter is due to confusion of sources with random polarization angles within a single beam element; 
although on average fainter sources with random (p will average to zero, their presence will introduce extra fluctuation into the 
polarization angle of the brightest source. 



A.4.2. Recovery of Source Counts 

Of particular interest to source surveys are the source counts dN/dS; it is therefore important to address whether this quantity 
can be accurately recovered in the presence of noise and a diffuse background. Figure 16 shows S^ '^dN/dS as a function of source 
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flux, for the input and recovered source distributions, in each of the simulation types Siml-Sim4. Since the input counts were 
dN/dS oc in this plot perfectly recovered counts appear as a line of zero gradient. The Figure shows that for Siml-Sim3, the 
counts obey the expected property of being well -recovered at high flux, but falling off as the flux approaches the survey detection 
threshold. Note that the fall-off at high flux at 150 GHz is due to the finite maximum 100 GHz source flux convolved with the 
assumed spectral index distribution; this part of the input counts does not obey a power-law but is well recovered by the source 
extraction. The source counts show similar properties at both frequencies, in total and polarized intensity. 
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Fig. 16. — Source counts S^'^dN/dS from simulated point source distibutions. Input source counts are gray, solid black is point-source only simulations (Siml 
and Sim3), dashed black is simulations including a diffuse background (Sim2 and Sim4). Top row is total intensity, bottom row is polarized intensity. Left pair 
of columns is 100 GHz, riglit pair is 150 GHz. Within each pair of columns, the left column is randomly distributed sources (Siml and Sim2), right is coirelated 
source locations (Sim3 and Sim4). The steepening of the 150 GHz counts at high flux is due to the convolution of a power-law flux distribution at 100 GHz with 
a gaussian spectral index distribution. 



Only in the case of the total intensity fluxes of clustered sources in the presence of a diffuse background (Sim4) do we see 
significant deviations from the ideal behavior; in this case, boosting occurs from low fluxes to higher fluxes. This may be seen 
from the defecit of sources at low flux relative to simulations without a diffuse component, and an excess above the input counts 
at high fluxes. In the high flux regime, where source counts are typically fit, the shape of S^'^dN/dS is heavily distorted from 
a simple power-law, rendering constraints on this quantity difficult to measure. Note the distinction between this type of flux 
boosting, which is due to background emission, and that due to un iform survey noise; in extragalactic surveys of radio sources, 
only the latter is normally considered (e.g. Muchovej et al. |2010] l, and can be corrected by marginalizing over the underlying 
source count parameters. For the former, which is termed 'background boosting' here, de-boosting source fluxes is not an easy 
problem. While diffuse emission such as the Cosmic Microwave Background (CMB) is very well characterized as a gaussian 
random noise distribution (and could potentially also be marginalized over to find the true source flux), the properties of the 
diffuse galactic emission (such as morphology, spectral behavior, and projection effects along different lines of sight through the 
galaxy) are at present poorly constrained and not easy to model. Contamination by galactic emission is also asymmetric in the 
sense that sources are only ever boosted to higher fluxes, because unlike CMB fluctuations (unpolarized) galactic emission is 
always positive. 

We therefore caution against over-interpretation of the source count slope in the presence of a galactic 'background', since the 
slope is shown by simulations to be corrupted despite the aggressive background filtering. The level of corruption is dependent 
on the parameters of the diffuse background model (such as the amount of power in beam-scale substructure). As we only use 
one set of parameters in the simulations, the results presented above are not intended to precisely quantify this systematic error, 
but to explore how the source counts can be affected by the background in a restrictive region of diffuse model parameter space. 
Despite these limitations, the source counts of the QUaD survey data shown in Figure |5] do appear to obey a power-law at high 
fluxes, indicating that in this regime background contamination is likely unimportant. 

A.4.3. Recovery of Spectral Index Distribution 

Figure[T7|shows the recovered total intensity spectral index probability distribution Pr{ai) for one realization of each simulation 
type. The input distribution, a gaussian p.d.f of unit rms, is well recovered for randomly distributed and spatially correlated 
sources without a diffuse background (Sims 1 and 3). Sims 2 and 4 exhibit a small shift of the distribution center towards larger 
values; this is likely due to contamination of faint source fluxes by the diffuse emission. In support of this notion, the shift is 
larger for Sim4 than Sim2 since more sources lie close to the plane of the galaxy, where diffuse emission is brightest. 

We conclude that the recovered spectral index distribution from the QUaD data (Section 4.4 1 may be biased slightly high as a 



result of background contamination preferentially affecting the 100 GHz data. However, the degree to which Pr{a) is corrupted 
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in simulations depends on the diffuse model input parameters; we therefore refrain from quantifying the effect. 
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AAA. Recovery of Correlation Function 

Simulations Sim3 and Sim4 are used to test the recovery of the input source spatial distribution parameters for sources clustered 
in the plane of the galaxy, as described in Section 4.5 The angular correlation function is constructed as in Equation]?] except this 
time Hci(6) (the number of sources with a neighbor at separation 6) is derived from a realization of Sim3 or Sim4, depending on 
whether a diffuse background is included. As with the real data, Siml is used to generate a histogram of the number of randomly 
distributed sources with a neighbor at separation 6, Hr{9). Correlation functions are generated for clustered point sources with 



and without a diffuse background present; the results for a single simulated realization are presented in Figure 18 
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Fig. 18. — Recovery of input source correlation function w(8) as a function of source separation 8. Input from one realization of the galaxy model is shown 
as black circles, with the recovered function shown as crosses for one realization of Sim3 (blue) and Sim4 (red), with small horizontal offsets applied for clarity. 
Left is 100 GHz, right is 150 GHz. The deviation from a power-law at 9 > 2° is due to the anisotropic source position distribution; this also occurs in the input 
simulations. At 6 < 0.4°, the power-law distribution is not well recovered — see text for an explanation. The recovered values of the power-law slope between 
0.4° <6 <2° are shown in the upper right of each panel. In this range, where survey shape effects are insignificant, the input slope is well-recovered. 



The black input points show that while a power-law correlation function is traced at small 9, at large source separations the slope 
becomes steeper This is a result of the anisotropic nature of the source distribution; less than 1 % of sources at each frequency lie 
further than 3° from the plane, and those inside are preferentially located towards b = Q (see Section 4.2 1. Only source clustering 
in galactic longitude contributes to the probability of finding a source separation greater than ^ 3°, suppressing the correlation 
function at large angular scales. The reconstructed w{9) from the simulated data shows this effect, demonstrating that although 
the correlation function is intrinsically suppressed due to the distribution of sources in the galaxy, it is still well-recovered at large 
angular separations. At angular scales < 0.4°, w(9) becomes poorly recovered. This is due to the large probability that the closest 
neighbors to bright, rare sources are faint and likely below the detection threshold. Therefore w{9) is only well-recovered when 
the source separation is large enough that the probability of a bright neighboring source is significant. 

Figure 18 shows that the slope of the correlation function is well-recovered in the range 0.4° < 6* < 2° to within the uncertain- 



ties; this is the range chosen for fitting w{9) in the QUaD data (Section 4.5 i 
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A.5. Effect of Background Kernel (Thck 

Without removal of the diffuse background, the measured flux of each extracted source in the QUaD survey can be heavily 
influenced by proximity to other bright sources and/or diffuse emission; either can add excess signal when source fluxes are 
determined, biasing recovered fluxes high. The filtering scheme described in Section [3] is designed to suppress the background 
by subtracting a template of diffuse emission from the maps. The template is a smoothed version of the raw survey map, with 
point source pixels replaced by their local median. Constructing the template requires a choice of smoothing scale abck, which 
represents the minimum angular scale on which background fluctuations are assumed significant. The results from the QUaD 
survey and the simulations in Appendices |A.2f|A!4| are dependent on the choice of (Ti^jt; here we investigate the effect of varying 
the value of this parameter. 
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Fig. 19. — Effect of o-jj-t on recovered total intensity flux S for simulated point sources above 25 Jy at 100 GHz. Solid lines coirespond to simulations without 
a diffuse background, while dashed lines include a model of diffuse emission. Left: Simulations of randomly distributed point sources. Right: Simulations 
of clustered point sources. For each simulation, the 16th, 50th and 84th percentiles of Sout/Si„ are plotted as a function of ctbckl '^heam- The solid vertical line 
corresponds to the value used on the QUaD data, cTbck/'^beam = 2.5, a compromise between excess filtering of source flux (low ai,ck} and excess background 
contamination (large abck)- 



Figure 19 shows the ratio of recovered 100 GHz / fluxes to the input as a function of (Jhck/ o'heam for sources above 25 Jy (similar 
results areround at 150 GHz). The left panel shows that for Siml, source fluxes are recovered to within ^ 2% if (Jtck/o'beam > 3, 
with a defecit of 5% by (Tbck/o'beam = 2. A similar result is found for Sim2, though increased variance is found for <Jbck / o'beam > 3. 
The systematic reduction of flux at low (Jbck/o'heam is due to source flux being subtracted with the background as ahck crheam- 
Increased variance is found when source locations are correlated, as in Sim3 (see right panel of Figure [T9|. Since this simulation 
is devoid of background emission, the larger 84th percentile is due to source confusion, as the likelihood of more than one 
source per beam is larger for spatially correlated sources. As might be expected, the 16th percentile is not significantly changed 
compared to randomly distributed sources since confusion cannot reduce recovered source fluxes. We note that despite source 
confusion, fluxes are still generally recovered to within 5% or better, independent of (Jbck- Including a diffuse background (Sim4) 
results in variance of S„ut/Si„ which increases as a function of cr bck / beam- The variance is skewed towards positive fluctuations 
due to the positive signal from the diffuse emission, which increasingly contributes to Sout as abck rises (because less diffuse 
emission is removed in the background subtraction stage). By (Jbck/ <^beam = 8, the 84th percentile of Soui/Si„ is ^ 1.2, compared 
to the 50th percentile of ~ 1.05; the fluctuation towards larger values of 5o,„/5,„ is therefore a+ = 0.15 or 15% higher than the 
input value. It is likely that this effect is also present in Sim2, where a diffuse component is present and sources are randomly 
distributed, but since far fewer sources are located close to the galactic plane, where the background is brighter, the effect is less 
obvious. 

The choice of <Jhck is therefore a compromise between loss of source flux due to excessive background subtraction (smaller 
abck), and increasing contamination from diffuse emission (larger abck)- Adopting abck/ abeam = 2.5 as in the QUaD survey data 
results in < 5% systematic loss of source flux, while reducing the scatter due to the background to ^ 5% — both these effects are 
comparable to the absolute calibration uncertainty in the maps of 3.5%. 
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